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We study the dynamics of a scalar field in the brane cosmology. We assume that a scalar field 
is confined in our 4-dimensional world. As for the potential of the scalar field, we discuss three 
typical models: (1) a power-law potential, (2) an inverse-power-law potential, and (3) an exponential 
potential. We show that the behavior of the scalar field is very different from a conventional 
cosmology when the energy density square term is dominated. 
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I. INTRODUCTION 

A scalar field has played a very important role in cos- 
mology. The problems of the standard Big-Bang theory, 
such as the horizon, and flatness problems could be re- 
solved by the inflationary scenario Q. In this scenario, 
the universe expands quasi-exponentially because of the 
vacuum energy, whose origin in many models are given 
by a scalar field (inflaton). Among inflationary scenarios, 
the chaotic inflation, which has a power-law type poten- 
tial (V{<j)) ~ 4> n : n is an integer) is the simplest model. 
Another example in which, a scalar field plays an impor- 
tant role is a cosmological phase transition and forma- 
tion of topological defects Q . The observed dark energy 
could be also explained by a scalar field whose energy 
density decreases in time (quintessence) ||[|. In a typ- 
ical quintessence model, the potential of the scalar field 
is assumed to be an inverse-power-law (V(<f>) ~ cf)~ a : 
[a > 0]). An exponential-type (V(<fi) ~ exp[— A0J) is nat- 
urally expected in many supergravity theories H . Hence 
those scalar fields are extensively studied by many au- 
thors in the conventional Einstein gravity theory. How- 
ever, recent higher dimensional unification scenario based 
on a brane world suggests that the gravitational law could 
be different from the Einstein's one at the early stage of 
the universe or at high-energy scale. 

In the brane-world scenario, which is based on a super- 
string or M-theory, our universe is embedded in higher di- 
mensions. Standard- mo del particles are confined in this 
four-dimensional hypersurfaces (3-branes), while grav- 
ity propagates in a higher-dimensional bulk space |6|-|9|. 
Among them, the Randall-Sundrum's second model is 
very interesting because it provides a new type of com- 
pactification of gravity. Assuming a 3-brane with a pos- 
itive tension is embedded in 5-dimensional anti de-Sitter 
bulk spacetime, the conventional four-dimensional grav- 
ity theory is recovered in low energy limit, even though 
the extra dimension is not compact U . While, at high en- 



ergy scale, gravity differs from the conventional Einstein 
theory |lC|j . Many authors discussed the geometrical as- 
pects of the brane and its dynamics (For a review, see, 
pl|), as well as its cosmological implications |l2] [T^]. 

Based on the Randall-Sundrum's second model, if we 
assume that our brane is homogeneous and isotropic 
(Friedmann- Robertson- Walker universe), the difference 
from conventional cosmology appears as two new terms 
in the Friedmann equation, i.e. the quadratic term of 
energy-momentum and dark radiation Jl2| , ^0| . As for the 
dark radiation, its effect is constrained by the nucleosyn- 
thesis. It is also diluted away in the inflationary era. 
Hence, the important change in scalar field dynamics, 
which we are going to discuss, may be due to the appear- 
ance of the term of quadratic energy density . 

The purpose of the present paper is to study dynam- 
ics of a scalar field confined on the brane, taking into 
account the effect of the quadratic-energy-density term. 
We discuss three typical potentials appeared in many cos- 
mological implications; a power-law potential, an inverse- 
power-law potential, and an exponential potential. The 
analysis of the quintessential potential will complete our 
previous discussion jl5|,|l6| . 

In Sec. [n[ we present the basic equations assuming the 
Randall-Sundrum's second model. Next, in Sec. EH, we 



show cosmological solutions in both scalar-field dominant 
and radiation dominant eras. These solutions are found 
analytically by using some approximations. 

In Sec. IV, we first study linear perturbations of the 
solutions, and then in Sec. we analyze those global 
stability. Sec. VI is devoted to summary and discussion. 



II. BASIC EQUATIONS 

We analyze dynamics of a scalar field in the Randall- 
Sundrum's second brane scenario Q , because the model 
is simple and concrete. It is, however, worthwhile not- 
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ing that the present result may be also valid in other 
types of brane world models, in which a quadratic term 
of energy-momentum tensor generically appears. In the 
brane world, all matter fields and forces except gravity 
are confined on the 3-brane in a 5-dimensional spacctimc. 
The extra-dimension is not compactified, but gravity is 
confined in the brane, resulting in the Newtonian gravity 
in our world. Since the gravity is confined in the brane, it 
is described by the intrinsic metric of the 4-dimensional 
brane spacetime. By use of Israel's junction condition 
and assuming Z 2 symmetry, the gravitational equations 
on the 3-brane is given by p0| 



(4) G„ = J4) Av + «4^ + 



E„ 



(2-1) 



where ' 4 'G Mt/ is the Einstein tensor with respect to the 
intrinsic metric <? M „, ( 4 ^A is the F 4-dimensional cosmo- 
logical constant, represents the energy- momentum 
tensor of matter fields confined on the brane and n M „ 
is defined by its quadratic form. E^ v is a part of the 
5-dimensional Weyl tensor and carries some information 
about a bulk geometry. k 4 2 = 8irG4 and k 5 2 = 8irGs are 
4-dimensional and 5-dimensional gravitational constants, 
respectively. In what follows, we use the 4-dimensional 
Planck mass m 4 = k 4 _1 = (2.4 x 10 18 GeV) and the 5- 



-2/3 



which could be 



dimensional Planck mass 7775 
much smaller than 771,4. 

Assuming the Friedmann- Robertson- Walker spacetime 
in our brane world, we find the effective Friedmann equa- 
tions from Eq.(|2.lD as 
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where a is a scale factor of the Universe, H = a /a is its 
Hubble parameter, k is a curvature constant, P and p 
are the total pressure and total energy density of matter 
fields, respectively. C is a constant, which term describes 
"dark" radiation coming from E^ u Jl2| ]. In what follows, 
we consider only the flat Friedmann model (k = 0) and 
assume that ^ 4 ^A vanishes for simplicity. 

As for matter fields on the brane, we consider a scalar 
field 4> as well as the conventional radiation fluid, i.e. 
p = p<j, + p T , where p^, and p r are the energy densities 
of scalar field <p, and of radiation fluid. We consider only 
the early stage of the universe, at which the quadratic 
term is dominant and then matter fluid can be ignored. 

Although a 5-dimensional scalar field living in the bulk 
fl9f may also appear in a brane world scenario, in this pa- 
per we only consider a 4-dimensional scalar field confined 
on the brane. The origin of such a scalar field might be 
found as a result of condensation of matter fields confined 
on the brane such as fermions. 

Since the energy and momentum of each field on the 
brane are conserved in the present model, we find the 



dynamical equation for a scalar field as a conventional 
one, i.e. 



• dV 

ZHcj) + — = 0, 



(2.4) 



where V is a potential of the scalar field. The energy 
density of the scalar field is 



P0 = + V(<j>). 
For the energy density of radiation fluid, we have 
p\ + 4H Pl = 0, 



(2.5) 



(2.6) 



which is integrated as p r cx a~ 4 . 

As the Universe expands, the energy density decreases. 
This means that the quadratic term was very important 
in the early stage of the Universe. Comparing two terms 
(the conventional energy density term and the quadratic 
one), we find that the quadratic term dominates when 



p> p c = 12m 5 6 /m 4 2 



(2.7) 



When the quadratic term is dominant, the expansion 
law of the Universe is modified. For example, the expan- 
sion law in radiation dominant era is a a i 1 / 4 instead 
of the conventional one t 1 ' 2 . Since we know the behav- 
ior of the scalar field in the linear term dominant stage, 
which is the conventional cosmological model, in this pa- 
per we study the behavior of the scalar field only in the 
quadratic term dominant stage. 

The nucleosynthesis gives a constraint on fundamen- 
tal constants, if the quadratic-energy-density dominant 
stage exists. The nucleosynthesis must take place in 
the conventional radiation dominant era to explain the 
amount of light elements. Hence, assuming that the 
energy density at a — a c is dominated by radiation 
as p c ~ (7r 2 /30)(?T 4 , where g is the degree of free- 
dom of particles, the temperature of the universe T c 
must be higher than that of nucleosynthesis, i.e. Tc > 
Tns ~ lMeV. This constraints implies 7775 > 1.6 x 
10 4 ( 5 /100) 1 / 6 (T NS /1 MeV) 2 / 3 GeV. 



III. COSMOLOGICAL SOLUTIONS 

There are two interesting limiting cases: one is the 
scalar-field dominant era and the other is the radiation 
dominant era. We discuss those solutions separately. 

A. Solutions in scalar-field dominant era 

Assuming the scalar field dominance, we find the Fried- 
mann equation (2.2) as 



H 



1 



6777 5 3 



-6 2 + V(<j>) 



(3.1) 
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We consider the following three types of the potential: 



V{4>) -- 

order. 



a+4, 



, M 4 exp[-A^] and ±m 2 2 , (|A0 4 ) in 



A.l V(4>)=n 



— ,, a +4j.-a 



The equation for the scalar field (2.4) is now 



3H — otfi' 



a+4. 



0. 



(3.2) 



Inserting Eq. (3.1) into Eq. (3.2), we find a second order 
differential equation for 0. It is easily shown that the 
asymptotic behaviors of the solutions are classified into 
three cases: (a) slow rolling (a potential term dominant) 
solution (a < 2), (b) a solution in which the potential 
term balances with the kinetic term (a = 2), and (c) a 
kinetic-term dominant solution (a > 2). 

Assuming a slow rolling condition, we find that 
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which lead to 



2am 5 3 . Then we find the solution 



= 2 [am 5 (t - t )] 

m 5 



1/2 



(3.4) 



where t is an integration constant. In order to keep the 
slow-rolling condition, we obtain a < 2 because (~ 

t-3/2) < t -( a +l)/2) and 02(_ t _l) ^ y-^ 

The universe expands as 



a = a cxp{[i/ (t-t )] (2 - Q)/2 }, 
where a constant i/o is given by 
Ho 



m 5 



= [3 (2 -a) (2^)"]^ 



JL 

m 5 



2(a+4) 
(2-e.) 



(3.5) 



(3.6) 



This solution describes an inflationary evolution, whose 
expansion rate is slower than that of the conventional 
exponential inflation, but faster than the power-law type. 
It may be interesting to discuss a spectrum of density 
perturbations for such an inflationary scenario. 
For a — 2, we find a power-law solution jTa] as 



a = a (t/t ) p , 



with 



P= 6 



l+~ 



JL 

m 5 



and 



m 5 



= 2[2m 5 t] 



1/2 



(3.7) 



(3.8) 



(3.9) 



The scalar field energy density evolves as 



oc t 



oc a 



-i/p 



(3.10) 



If p > 1, i.e. [i > 40 1 / 6 m5 ps 1.85ms, we have a power- 
law inflationary solution |20|] , which is an attractor of the 
present system as we will see later. While, if p < 1/4, 
i.e. \i < 4 1 / 6 m5 w 1.26m5, this solution is no longer an 
attractor, leading to the radiation dominant era (see Sec. 

If a > 2, a kinetic-term dominant solution gives the 
asymptotic behavior of the present system (see Sec. IV). 
If a kinetic term is dominant, we find 



H 



leading to 



m 5 
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and + ZH<t> « 0, 



(3.11) 



(l/4m 5 3 )0 3 , which is integrated as 



±2v / 2m5 /2 (t - t ) 
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(3.12) 



where to an d 0o are integration constants. For the 
solution with + sign, the potential term decreases as 
V oc t~"/ 2 , while the kinetic term drops as 2 oc t -1 , 
then the kinetic term is always dominant if a > 2. In 
this case, we find that the Universe expands as 



a = aa(t/t ) 1/6 , 



(3.13) 



which is exactly the expansion law of the quadratic-term 
dominant universe with stiff matter. 

On the other hand, for the solution with — sign, 0^0 
as t — > to 5 and then V — > oo. Then before reaching to 
, the potential term becomes dominant. However, since 
the potential-dominant solution is not any attractor, it 
turns out that the former kinetic-dominant solution is 
eventually reached (see Sec. IV). Consequently, the 



kinetic-term dominance is the asymptotic behavior for 
a > 2. This guarantees the radiation, if it exists, eventu- 
ally dominates the scalar-field, which makes initial con- 
ditions for a successful quintessence wider [fMTF 



A.2 7(0) = ^ 4 exp[- 



■A-2-1 

7715 J 



The exponential potential is steeper than the inverse- 
power-law pot entia l. Th en the asymptotic solution is the 
same as Eqs. ( 3.12 ) and ( 3.13 ) in the case (1) with a > 2. 
Then the radiation dominant era is eventually realized. 

A. 3 Chaotic inflationary potential 

The inflationary solution of this model was already ob- 
tained (22), for completeness, we list here the analytic 
solution including its oscillating phase. 

For power-law potential such as V — \m 2 <\? or |A0 4 , 
we find a stronger inflation than the conventional one 

^For the model V(0) = i?n 2 2 , assuming a slow rolling 
condition, we find that 
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12ml 



and 



m 
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(3.14) 
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which lead to 



-4m 5 3 . We then obtain the solution 



±2v / 2m 5 3/2 (*o-i) 1/2 , 



(3.15) 



a oc t 1/4 and (</>) cx t~ 1/4 . 



(3.22) 



Introducing a conformal time r\ = J dta(t) 1 and the 
conformal field tp = acf>, the equation for tp (Eq. (EES)) is 



a = ao exp 



(3.16) 



where to and ao are integration constants. 

The slow-rolling condition is broken when \(f>\ ~ |V'| 
or (l/2)0 2 ~ V. Both conditions give the time when 
inflation ends, i.e. tf ~ to - (2m) . In order for inflation 
to take place in the quadratic energy density dominant 
stage, P4,{4>f) > p c , which implies m > im^/m 2 . 

After the end of inflation, there appears an oscillation 
period, which solution is given approximately as 



2m 



3/2 



mi 1 / 2 



sin mt, 



a = a (t/to) 1/3 . 



(3.17) 



(3.18) 



The scale factor evolves just as the case with a pres- 
sureless perfect fluid (dust fluid) at the quadratic term 
dominant stage. Since the amplitude of <p decreases more 
slowly than the conventional case because of the small ex- 
pansion rate, we expect a sufficient particle production 
in the preheating phase Q . 

(ii) For the model with V(<p) = jXcj) 4 , the behavior of 
the scalar field in the inflationary phase is not so different 
from the case of the previous massive inflaton model. 
Assuming a slow rolling condition, we find that 
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which lead to 
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a = ao exp 



9 



(to - tf 



(3.20) 



(3.21) 



where to an d «o are integration constants. 

The slow- rolling condition is broken at t — tf ~ 
t - 2 1 / 3 /(A 1 / 3 4m 5 ) when (l/2)0 2 ~ V. Another con- 



dition 



|V'|) corresponds to the time tj 



to - 



l/(A 1 / 3 4rri5), which appears after tf. The scale factor at 
the end of inflation cif = a(tf) — aoexp[— 1/9] ~ ao- 
In order inflation to take place in the quadratic en- 
ergy density dominant stage, p<f,((j)f) > p c , which implies 
A > (27/32) (m 5 /m 4 ) 6 . 

For the oscillating period, it is shown that there is a 
time-averaged relation between the potential term and 
the kinetic term of the scalar field, i.e. (</> 2 ) = A(</> 4 ) 
and then the scalar field behaves approximately as the 
radiation fluid Jl7| . Then we find that 



ip" + \<p 3 tp = 0, 

a 



(3.23) 



where ' stands for the derivative with respect to the con- 
formal time. Since the time averaged value of ip is almost 
constant while (a" /a) ~ f] 2 ^ the third term can be ig- 
nored. In the conventional cosmology, since a ~ t 1 / 2 ~ 77, 
the third term exactly vanishes. Then we can integrate 
Eq.(3.23), finding the elliptic function as 



tp = <p en l x — xq 



V2 



(3.24) 



where a constant <p is an oscillation amplitude of tp, x = 
^/\tprj is a dimensionless conformal time variable, and Xq 
is the value of x when the inflaton field starts to oscillate 
coherently. This solution is the same as that obtained in 
the conventional cosmology. However the corresponding 
time evolution of the time-averaged scalar field (<fr) is 
given as {<fi) cx a -1 oc t -1 / 4 , which decreases more slowly 
than the conventional one. We study this effect at the 
preheating stage in |lq|. 



B. Solutions in radiation dominant era 

In the radiation dominant era, the scale factor expands 
as a oc t 1 / 4 . We have to solve the equation for a scalar 
field. 

The basic equation is now 



(3.25) 



We shall consider again the same potentials discussed in 
the previous section in order. 



B.l V(^)=fi 



..a+A^—a 



The equation of motion for the scalar field (2.4) is now 



3 ^ 
— © — au 

4t 



a+Ai—a—l 



= 0. 



We find an analytic solution for a < 6 [ Jl5| |, that is 

2 

t \ °+ 2 



with 



i«+2 



2a(a 



6 — a 



-n a+A t 2 



(3.26) 

(3.27) 
(3.28) 



where to and 4>o(> 0) are integration constants. Eq. 
Q3.28 ) requires a < 6. 

The energy density of the scalar field evolves as 



4 



3a(a + 2) TA / t 
P4> = Vo — 

3a(a + 2) / a 

- — ~ r Vo — 

D — a \ao 

where Vo = V(0o) and ao = a(to). 

The density parameter of the scalar field, which we denote as fi^, is 



(3.29) 



Pr + P4, Pr 



«0 



(3.30) 



where fi^ = 3 "g 2 ^ Vo / p r (t Q ) is an initial value of the density parameter. 

If a > 2, fi^ decreases with time, just contrary to the tracking solution in the conventional quintessence model. The 
scalar field energy decreases faster than that of the radiation. This is a new interesting feature because the "initial" 
smallness of a quintessence-field energy could be dynamically explained If a — 2, is constant until the 

liner term becomes dominant. This is the so-called "scaling" solution. The scalar field energy drops at the same rate 
as that of the radiation. On the other hand, if a < 2, increases, and then the universe get into a scalar field energy 
dominant era discussed above, which is an inflationary solution. 
For a > 6, we find a kinetic-term dominant solution 



1/4 



(3.31) 



where </>o(> 0), and to are integration constants, decided by the initial value only. In this case, they do not include the 
characteristic scale, because there is no contribution of the potential term for the energy density of the scalar field. 

B.2 V(0)=/i 4 exp[-A^-] 



This case is not so interesting, because the kinetic energy dominant solution ( |3.31 ) gives an asymptotic behavior 
just as the same as the case of scalar field dominant era. 

B.3 Chaotic inflationary potential 

(i) For the model V((f>) = ^m 2 (j) 2 , the solution is given approximately as 



>J0 \r 



-3/8 



sin rat. 



(3.32) 



The energy density of the scalar field p^, which evolves approximately as cx t~ 3 / 4 cx a~ 3 , is the same as that of 
dust fluid, then the scalar field energy will eventually dominate unless its energy is transferred into radiation via a 
reheating process. 



(ii) For the model V{4>) — jX(p 4 , the behavior is the same as Eq.(3.24), because the time evolution of the scale 
factor is the same. 

We summarize the above solutions in the next table. 



R (a cx t 1 / 4 ) 



potential 



a < 2 



a = 2 



2 < a < 6 



a > 6 



k 2 ^ 



inflation : Eqs.(3.5|),(|3.4 



power-law : Eqs . (|3 . 7|) , (|3 S§ 



kinetic dominance: Eqs.( |3.13 ),( 3.12D 



power-law : Eq.(3.27) 



kinetic dominance:Eq.(3.31) 



oscillation:Eq.(|3.32|) 



inflation:Eqs.(3.1£),(3.15 



oscillation: Eqs . ( 3 . 1 



(3.17) 



inflation:Eqs.(3.21),(3.20) -> oscillatiomEqs. (3.22), (3.24) 



oscillation:Eq.(3.24) 



TABLE I. The asymptotic behavior in quadratic-energy-density dominated stage. S and R denote scalar-field dominant era (p^> 2> pr) 
and radiation dominant era (p r 2> p</>), respectively. The arrows describe the final fate. 



5 



IV. STABILITY ANALYSIS 

In the previous section, we have listed up several ana- 
lytic solutions in the quadratic energy density dominated 
era, which are summarized in Table 1. If those solutions 
are one of the attractors in the dynamical system, such 
spacetimes may be naturally realized in the history of the 
universe. Therefore, it is very important to study their 
stability. 

In this section, we study the stability of the solutions 
against linear perturbations. Although we may need the 
detailed analysis for the solutions for the chaotic infla- 
tion and the following oscillation, we expect that such 
solutions are found for a wide range of initial conditions 
just as the same as conventional chaotic inflationary sce- 
nario. Hence we focus only on the stability of solutions 
of inverse-power-law and exponential potentials. 

Since the cosmological solutions are time dependent, 
we first have to find new variables by which the solutions 
are described as fixed points in the dynamical system. 



A. scalar field dominated stage 

1. the power-law solution 



First, we investigate the power-law solution (3.7), (3.S) 
with an inverse-power-law potential with a = 2, i.e. V = 
/i 6 0~ 2 . The asymptotic behavior of the solution is 



H 

771 5 

V 

4 

mt 



P 



(m 5 t) 
<o 2 



(4.1) 
(4.2) 
(4.3) 



where M = (fi/msf /2y/2 and p = (1 + M 2 )/6. 

For these time dependence, we introduce new variables 

as 



H = tH, 
V = mJ 3 tV, 



(4.4) 
(4.5) 
(4.6) 



In terms of these variables, the basic equations are ex- 
pressed as, 



W =H(1-6H + V), 

V = V - (V27Wo) _1 PV 3/2 , 

V = -(1 - m) + (V27W ) _1 V 3/2 , 



with the following constraint, 



n = T2 V 



iv. 



(4.7) 
(4.8) 

(4.9) 
(4.10) 



where a prime denotes a derivative with respect to new 
time coordinate r = ln(77i5i). 

These suggest that there are only two independent 
variables. 

For the basic equations, the power-law solution is de- 
scribed by a fixed point (Hp, Vq, Vq) — (j>, M 2 

After eliminating TC with Eq.(4.10), the perturbation 
equations around this solution are 



SV 



Ml 
V2 



5V; 



SV. 



(4.11) 
(4.12) 



Because these equations do not include the explicit 
time-dependence, we can discuss the stability of the so- 
lution, by the eigenvalues of these coupled equations. 

They are 



(4.13) 



which are all negative, and then this solution is stable 
against linear perturbations. 



2. the kinetic dominant solution 

Next, we analyze the stability of the kinetic dominant 
solution ( p. 12 ) with an inverse-power-law potential with 
a > 2 and an exponential potential. First we discuss the 
case with the potential V(<fi) — fi a+4: <fi~ a . The asymp- 
totic behavior of the solution is 



H 

m 5 
V 

4 

mi 



6(m 5 t) 

V Q (m 5 ty^ 2 , 

A. ^ 

m\ ~ (m 5 1/2 ' 



(4.14) 
(4.15) 
(4.16) 



where V = 8- q / 2 (/j/to 5 ) q + 4 . 

For such time-dependenc e, th e variables TL and V are 
given by the same relations (4.4) and (4.6), but V is newly 
introduced as 

V = TOjTVe^ r , (4.17) 

using the time coordinate r = \T\{m^i). Then we find the 
following basic equations as, 



W = H{l-m + Ve-^ T ), 



V 



-V 



V 



p' = _(l-6W) 



V 

a 



2V2 ° 



(4.18) 
(4.19) 

(4.20) 
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with the following constraint equation, 



H 



1 



Ve 



(4.21) 



From Eqs. fl4.18|) - fl4.20D with ( |4.2l| ), in the limit of 
t — > oo j we find the fixed point (H, V,V) = (1/6, V , \/2), 
which corresponds to the kinetic-term dominant solution. 
Eliminating Ti with Eq.(4.21) and expanding the above 
basic equations around this fixed point as {7i,V ,V) = 
(1/6 + SH, V + SV, V2 + SV), we find that the perturba- 
tion equations are 



2 2\/2 

SV' = {a- 2)-^e~^ r - SV. 
2v 2 



(4.22) 
(4.23) 



One may think that this system is easily analyzed be- 
cause the time-dependent part in these equations decays 
with time. If we estimate the eigenvalues for the asymp- 
totic equations, we find that those are -1/2 and -1, which 
suggests the system is stable. However, since the equa- 
tions have explicit time-dependence, we have to be care- 
ful to conclude it Fortunately, in this case, we can easily 
integrate the above perturbed equations directly as 



SV = 5V e~% + ^SV a e~ T 
v 2 



a(«-2) _ v , ^ 



2(a-3)(a-4) 
SV = SV e- r - t~ 2 \, Voe-' ~ 



y/2(a - 4) 



(4.24) 
(4.25) 



if a ^ 3,4, 



3V 

vr 



SV = ^SV e^ T - -V 2 re-2 T + SV e~^ T (4.26) 



SV = 5V e- T + — V e-* T , 
for a = 3, and 



(4.27) 



SV = 8V Q e-^ T + 2V2V SV Q e- T + 2V 2 (t + 2)e~ T (4.28) 

SV = SV e- T + -^-V Q re- T (4.29) 
v2 

for a = 4, where SVq and SVq are the integration con- 
stants. These show that SV, SV — ► in the limit of r — > 
00, which means that the kinetic dominant solution is 
stable against linear perturbations for V((j)) — fi i+a (f>~ a , 
with a > 2. Note that for the case of 2 < a < 3, the 
perturbations drop more slowly than those expected just 
from the above eigenvalues. 

Next, we consider the case with an exponential poten- 
tial V = (J, 4 e~ X( t '• tm '~> . The asymptotic behavior of the 
solution is 



1 



H 

V2 



^-^2V2(m 5 t) 1/2 , 
m 5 



to 2 (m^t) 1 / 2 



(4.30) 
(4.31) 
(4.32) 



In this case, because the potential decreases quite fast as 
V oc exp(— 2v / 2A(m5t) 1 / 2 ), we adopt new variables 



T 



(4.33) 



instead of V. The basic equations are given as 

4 



H' = H 



1 

24' 



4 



JL JL 

12 \m 5 



V 2 exp (t - Xe r/2 T 



T' = --T 
2 



V 



|(1 - QH) + X^ exp Qr - Xe T ' 2 T 



(4.34) 
(4.35) 
(4.36) 



with the following constraint equation, 



(4.37) 



If T > 0, in the limit of t — * 00, we find a fixed point 
{Ho,Tq,Vq) = (1/6, 2\/2, V%), which c orres ponds to the 
kinetic-term dominant solution ( |3.12j ),( pTl~3 ). 

Expanding the basic equations around this fixed point 
as (H,T,V) = (1/6 + SH, 2V2 + ST, V2 + SV), and elim- 
inating H with Eq. (|4.37|) , we find that the perturbation 
equations are 



ST' 



1 



-ST + SV, 



SV = -SV + X^-r exp ( - 2 v / 2Ae r/2 



mi 



V2 " 4 



— j exp r 

2 mt 



2V2Ae T/: 



(4.38) 



(4.39) 



In this case the time-dependent term in the equations 
drops very fast as exp(— 2v / 2Aexp(r/2)), we just analyze 
the asymptotic equations In the limit of r — * 00, which 
is 



ST' = --ST- 
2 



SV, 



SV' 



-SV 



(4.40) 
(4.41) 



The eigenvalues of the system ( |4.40| ) , ( 141 ) are -1/2 and 
-1, then the kinetic dominant solution is stable against 
linear perturbations for V(<fi) — ^j 4 e~ A */ m5 as well. We 
can easily check that it is correct even including the time- 
dependent terms. 
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3. the inflationary solution 

Final ly, w e study a stability about the inflationary so- 
lution (3.4) with the inverse-power-law potential V — 
fi a+4 (j)~ a (a < 2). We use H, V(4>) and <j> as our dynami- 
cal variables. The asymptotic behavior of these variables 
for the inflationary solution is 



H 

m 5 



V 

6m| 



\v Q (m 5 ty 

6 



Vo {m 5 t)~ 



(4.42) 
(4.43) 



where Vo = {2y/a)~ a (fi/m 5 ) a and Vo = a? . For such 
time-dependence of the solution, we introduce new vari- 
ables as 



U 
V 



m - l H TP^T, 



V = ml 2 ^ T~ 



(4.44) 
(4.45) 
(4.46) 



where the new time coordinate T is defined as 

T = {m^t) 1 ^ . (4.47) 
Using these variables, the basic equations are now, 



If we adopt a naive analysis for the asymptotic equa- 
tions ignoring decaying time-dependent terms, we find 
that the eigen values are — Vo/(2 — a) and 0. The zero 
eigenvalue corresponds to a marginally unstable mode, 
which appears in a slow-rolling phase of an inflationary 
solution. 

In order to confirm the asymptotic behaviors of <5V and 
SV, we will integrate the above equations directly. To 
solve (4.52) and (4.53), we first need to find the following 
first-order differential equation 



— [A(T) SV + B(T) SV} = 

a \A(T) SV + B(T) SV] + F(T). 



(4.54) 



where A(T), B(T), and F(T) are chosen appropriately 
as follows. 

We can easily show that if A(T), B(T), F(T) satisfy 



dA 
dT 
dB_ 
~dT 



Vn 



A 



B 



V 



(2-a)a 



A 



(2 - a)T T 
1 \ B 



0, 



and 



F = 



V 



(2 -a)T J T 



1 2 \ A 
1 ~ 2 / T> 



0. 



(4.55) 
(4.56) 

(4.57) 



/ a\ dH 
V 2/ ~dT 



a\ dV 
~ 2 ) dT 



a\ dV 
2/ dT 



-H(<oH-V) + -H t-\ 



a 

2 V 



Vo \Vo 



l/a' 



vv m (v 

V -T-\ 
2 



with the following constraint equation, 



H 



1 

12' 



-V T~ 



-V. 



(4.48) 



(4.49) 



l/a 



(4.50) 



(4.51) 



( |4.54| ) is found . 

From ( 4.55| ) and ( 4.56] ), by eliminating A, we find the 
second order differential equation for B, as, 



d 2 B 
'dT 1 ' 



2a 



Vo 

2- a 

Vo 



1 1 



2 - a 



1 



(2-a)T 2 



x B 



2-aT 
a 

2~a + 
0. 



dB 

2V 
(2 -ay 



1 

T 

(4.58) 



The solution of this equation is given by two indepen- 
dent hypergeometric functions, which asymptotic behav- 
iors are given by 



In the limit of T — > oo, we find the fixed point 
(Vo/6, Vo,Vo), which corresponds to the inflationary so- 
lution. 

Expanding the basic equations around this fixed point 
as (H, V, V) = (Ho + SH, V + SV, V + SV), and elimi- 
nating H with Eq.(4.51), we find that the perturbation 
equations are 



1 - 
1 - 



a\ dSV 
2 J ~dT 
a\ dSV 
2J ~dT 



-\5VT- 1 - 

^sv-^ 

2a 2 



aVo 
2V 

SV 



5VT- 



1 , 1 9s 1 

2^0(1 - ^Pl)T-\ 



(4.52) 



(4.53) 



B 1 (T) 
B 2 (T) 



„rrTrp^- 
e 1 2 ~ c 



[<*+ 



v 



J 1 - — 



(4.59) 
(4.60) 



From (|4.56|) and ([4.571) , we find A(T) and F(T) for 
each solution as 



Fi(T) 
F 2 (T) 



2- a 
aVo k+ . 

Vo 

aVo 
2-a 



„ 2(Q-1) 



v 



Vn 



\T T - 



(4.61) 
(4.62) 
(4.63) 
(4.64) 
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We then obtain two set of (4.54), which general solu- 
tion is given by 

A X {T) SV + Bx{T) SV = e aT J dT'e^ aT ' Fi(T'), (4.65) 
A 2 {T) SV + B 2 {T) SV = e aT J dT' e - aT F 2 (T'). (4.66) 
The asymptotic behaviors of these equations are 

Q 



V 



aVo 

' v 



SV + T5V 



Cl T-v=z ( 1 - -VZ 



1 



(4.67) 



SV + 5V 



= c 2 T^e-^ T -a(l- 1 -V^T-\ (4.68) 

where c\ and c 2 are integration constants. 

Then we obtain the asymptotic solution of SV and SV 

as 



SV = &T l [ciT-^ - c 2 T^e 
+ (*-l)(l-\vi\], 



(4.69) 



SV 



. , . aVn „ a a i -vp rr 

A- 1 ci — -T _ — -c 2 — e — T 
1 7>o V 

a(a — VqT f 1 



V T 



l--Vr 



J7 



where 



a 



A=-[V T-1] 



(4.70) 



(4.71) 



( ^69| ), jl7C|) and (|jT]) give the asymptotic behavior 
of SV and <5V as 



SV 
SV 



(a - l)Vo 
aVn 



1 - -zVl T- 1 



(4.72) 
(4.73) 



This power-law decay corresponds to zero eigenvalue. 
Then the above naive analysis by eigen values is con- 
firmed. We conclude that the inflationary solution is sta- 
ble against linear perturbations. 



B. radiation dominated stage 

1. the power-law solution 

We show here that for the scalar field with an inverse- 
power-law potential model V(</>) = /x 4+Q 0~ Q ~ 4 , with 
a < 6, the power-law solution in the radiation dominant 



stage given by ( 3.27) , (3.28) is stable against the linear 
perturbations. 

Following 0], we introduce the new variables as 



r : 



(4.74) 
(4.75) 
(4.76) 



where <f> e is the exact solution given by ( 3.27 ) , ( 3.28 ) , 

— = A(a)(m 5 i)^T2. (4.77) 
m 5 

A(a) is a dimensionless constant, which depends on a 
and /X, 



A(a) 



2a(a + 2) 2 



- 

\m 5 



. I 2 



6 — a 

With these changes, Eq.( [3.26j ) becomes, 



T' = V, 

V = 



14 



a — 6 



4(a + 2) 2(a + 2) : 



a — 6 
2(a + 2) 



-T 



(4.78) 



(4.79) 
■ a_1 , (4.80) 



where prime denotes the derivativ e wit h resp ect to r and 
the power-law solution given by ( 3.27 ) and ( 3.27 ) corre- 
sponds to a fixed point (J 7 , V) = (1,0). 

In order to study the stability near the critical point, 
we here use the linear analysis Expanding the 

above equations around the fixed point as {T,V) — 
{l + 5F,5V), we find that 



ST' = SV, 

2(a + 2) 



a - 14 
4(a + 2) 



SV. 



(4.81) 
(4.82) 



Because these equations do not include the explicit 
time-dependence, we can discuss the stability of the so- 
lution by the eigenvalues of these coupled equations. 

The eigenvalues are given as 



and — 



6 — a 



(4.83) 



which are all negative for a < 6, and then this solution 
is stable against linear perturbations. 



2. the kinetic dominant solution 



We show here that for the scalar field with an inverse- 



4+ct J. — ol — 4 



power-law potential model V(<p) = fi 
6, and an exponential potential model V(<p) — fi 



with a > 

4 -\<t>/m 5 



the kinetic dominant solution in the radiation dominant 



stage given by (3.31), is stable against the linear pertur- 
bations. 







At first, we show the case with an inverse-power-law 



potential model V(<f>) = fi 



4+Q . 



, with a > 6. As in 



the power-law solution in the radiation dominant stage, 
we introduce the new variables as 



t : 
T - 



m 5 t, 



(4.84) 
(4.85) 



where cf>k is the kinetic dominant solution given by (3.31) 



m 5 



(4.86) 



In this case, B is the dimensionless constant, whose 
value depends on the initial c ondit ion. In terms of these 
variables, the basic equation ( 3.26 ) is expressed as, 



T' = V, 



m 5 



a+4 



(4.87) 



In the limit of r — > oo, we find the critical point 
(T, V) = {Tq, 0), where To is an arbitrary constant. This 
point corresponds to the kinetic dominant solution. Ex- 
panding the basic equations around this critical point as 
(T,V) = (Tq + ST,SV), we find that the perturbation 
equations are, 



ST' = SV, 



SV' 



1 



-SV 



m 5 



a+4 



zr 



a — 2 -j^—a — 1 
J~ n 



e 4 



(4.89) 
(4.90) 



Because 



these equations 
have explicit time-dependence, we cannot evaluate the 
stability with the eigenvalues. However, for a > 6, we 
can easily integrate the above perturbed equations as 



ST 



16a ( [i 
(7- a)(6 - a) \m 5 
A5V e-i + ST , 



a+4 



B- 



a — 2 sr-— a — 1 
J~ n C 



SV = 



4a 



B 



-a-2 rr — ot — 1 , 



7 — a \ ms 
+ 5V e~-- 

for a ^ 7, and 



<)T = -2S [ JlY 



(4.91) 



(4.92) 



S- 9 ^ - 8 e-*(r-4) 



4#> e"~* + ST , 
n 



SV = 7 



JL 

7715 



9 t- 8„ 



(4.93) 
(4.94) 



for a = 7, where & n d #Po are integration constants. 

These show that in the limit of r — > oo, even though 
this is not strictly a fixed point, the kinetic-term domi- 
nant solutions are attractors in a sense that these satisfy 

(Plb- 



Next, we investigate the case with an exponential po- 
tential model V((f>) = ^e"^^ 5 . 

In this case, the equation of motion for the scalar field 
( |2.4| ) is now, 

" >! ' ' ! " ' (4.95) 



At 



Xfj, — exp[— Xtfi/ms] 

7775 



0. 



In terms of the variables, ( 4.84 )- ( |4.86| ), the basic equa- 
tion (4.95) is expressed as, 

T' = V, (4.96) 



V 



V, 

1 



exp 



t - \Be? T T 



(4.97) 



A / n 
—V H — 

4 B \"i5. 

In the limit of r — » oo, we find the critical point 
(T, V) = (To, 0), where Tq is an arbitrary constant. This 
point corresponds to the kinetic dominant solution. Ex- 
panding the basic equations around this critical point as 
(T,V) = (Tq + ST,SV), we find that the perturbation 
equations are, 

ST' = SV, 

1 



SV 



-SV 



B 



JL 

777 5 



exp 



t - \Be* T T 



(4.98) 
(4.99) 



In the limit of r 
these equations are, 

ST' = SV 



oo, the asymptotic behavior of 



SV 



1 



-6V, 



r/4 



(4.100) 
(4.101) 

(4.102) 
(4.103) 



which can be integrated as 

ST = 5T - A5V Q e 

SV « SV e- T /\ 

where STq and SVo are integration constants. 

These show that in the limit of r — > oo, even though 
this is not strictly a fixed point, the kinetic term domi- 
nant solutions are attractors in a sense that these satisfy 

(fEU- 



V. GLOBAL STABILITY OF THE SOLUTIONS 

In the previous section, we show that the analytic so- 
lutions listed up in Table 1, are stable against the linear 
perturbations. However, in order to show that these so- 
lutions are realized from wide initial conditions, we have 
to study stability against non-linear perturbations. The 
we analyze a global stability of the solutions in the phase- 
space. 

For the model with the potential V((f>) — p!°$r 2 , we 
study it both numerically and analytically including both 
the scalar field and the radiation (§VA). We also show 
the numerical results for the kinetic dominant solutions 
( |3.12[ ) and the inflation ary s olutions (3.4) by assuming 
scalar field dominance (§VB) and for the power-law so- 
lutions (a 2) (3.27) and the kinetic domi nant solutions 
( |3.3l| ) by assuming radiation dominance (§ VC ). 
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A. the power-law solution (a = 2) 

Here, we study a global stability of the solutions for 
V(4>) = H 6 4>~ 2 in detail, following p4| . We introduce 
new variables 



Y' <o 



x= m 5 



3/2 1 



2V3 Vff' 
2ml' 2 r 1 

Y —jr Mo 7M' 



(5.1) 



where Air 



7(2v^r, 



and new derivative with re- 



spect to In a, which is described by '. The dynamical 
equations are rewritten as a plane-autonomous system: 

X' = F(X, Y) = X(X 2 - 2Y 2 - 1) + 37W _1 y 3 , 

Y' = G(X, Y) = -Y (2Y 2 - X 2 + ZM^XY - 2) , (5.2) 

and radiation energy is given by the constraint equation 



Pr 



X 1 



Y z = 1. 



(5.3) 



From the constraint equation, the density parameter 
of a scalar field is given by 



QmiH 



X 2 + Y 2 



(5.4) 



Since radiation energy is non- negative (p r > 0), we 
find < X 2 + Y 2 < 1, so the evolution of this system 
is completely described by trajectories within the unit 
disc. As the system is symmetric under the reflection 
(X, Y) — > (X, —Y) and the evolution will not go beyond 
the Y = 0-line (which corresponds to <p = oo or H = oo), 
it is enough to discuss the upper half-disc (Y > 0). 

Depending on the value of Moiorp/mc,), we have four 
or five fixed points (critical points) where X' and Y' van- 
ish. The four critical points (0,0), (1,0), (—1,0), and 
A(l/ v /l + M 2 , + M 2 ) exist on the bound- 

ary of the phase space, while the fifth critical point B 
(</870Mo /2 , {/2/9Ml /2 ) appears in the half unit disc 
only if Mo < V2/2. 

The simple analysis by linear perturbations shows that 
the critical point A is stable if A4q > v / 2/2, while it 
becomes a saddle point when Mo < \/2/2, and newly 
appeared critical point B becomes stable. Since the point 
A locates on the boundary, is always unity (scalar- 
field dominated stage), while the point B appears on the 
line Y = X/^/2. The density parameter of this solution 
is flff, = \/2Mo, which means that the radiation density 
becomes larger that that of a scalar field if Mo < v2/4 

In order to analyze the global behavior in the phase 
space of the present dynamical system, we first draw crit- 
ical curves corresponding to F(X.Y) — 0, which corre- 
sponds to the curve C\ or G(X, Y) = 0, which gives two 
curves;the straight line Y = and hyperbolic curve C2 
(see Fig. |and Fig.|). 




Magnified figure for the rigion around the point A 
Y 

0.1 0.2 0.3 0.4 

1.0 | ^^ ' ■ ■ 7 X 




FIG. 1. Schematic phase diagram for the model with 
V{4>) — ^ 6 <t>~ 2 for the case of Mo > The arrows show 

the directions in which the spacetime evolves. The critical 
point A is an attractor. 




Magnified figure for the rigion around the point A 




0.5 0.6 0.7 0.8 0.9 



2. Schematic phase diagram for the model with 

6j-2 



FIG 

V(4>) — Li 6 (j>~ 2 for the case of Mo < v/2/2. The arrows show 
the directions in which the spacetime evolves. The critical 
point A is a saddle point, while the point B is an attractor. 

From the signs of F(X.Y) and G(X, Y), We find the 
direction of the evolutionary track of the universe in the 
phase space, which is shown by arrows in Fig. [J and 
Fig.||. Following those arrows, we find whether the crit- 
ical points are globally stable or unstable. It turns out 
that the critical points (1,0) and (—1,0) are unstable, 
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and (0, 0) is a saddle point. For Mo > y/2/2, the critical 
point A is stable, but if Mo < \/2/2, it becomes a saddle 
point, as we expect from our perturbation analysis, and 
new critical point B becomes a stable point. 

We also show the result of numerical analysis for both 
cases. We set \i = 2m^ (Mq = 2^/2) for the case of 
Mo > V2/2 and /i = m 5 (Mo = \/2/4) for the case of 
Mo < V2/2. 



responding to the critical point B becomes an attractor. 



B. scalar field dominated stage 

In what follows, we only show numerical analysis for 
global stability for other cases. 



1.0 

>- 
0.5 



(7) 
0.0 
-1.0 



/ / 



if / + 

(6) (5)-,' (4f 
-0.5 0.0 

n/m 5 = 2 



\ s 

\\ I 

\ 

\ \- 

(3) (2) j (D 
0.5 1.0 X 



V. 



FIG. 3. Numerical analysis for the model V(4>) = fj 6 ^ 2 
with fx = 2m 5 (Mo = 2y/2 > %/2/2). The critical point 
(X,Y) = (1/3,2^2/3) corresponding to the power-law solu- 
tion (3.7) and (3.9) is the attractor. As for initial values of 
(</>/m.5, (4>/m$)' = t x d(4>/ms)/dt,pr/m\) at m^t — 1.0, we 
set (1) (1.0 xlO 2 , 1.0, 1.0 xlO -4 ), (2) (1.0 xlO 2 , 1.0, 1.0 
xlO -1 ), (3) (1.0 xlO 2 , 1.0, 1.0 ), (4) (1.0 xlO 2 , 0.0, 1.0), (5) 
(1.0 xlO 2 , -1.0, 1.0), (6) (1.0 xlO 2 , -1.0 , 1.0 xlO -1 ), (7) 
(1.0 xlO 2 , -1.0 , 1.0 xlO -4 ). 




A (6)/ 



0.0 

-1.0 



(i) 



-0.5 v 0.0 
/ u/m 5 = 1 



0.5 



1.0 



X 



FIG. 4. Numerical analysis for the model V(cf>) = p 6 cj> 2 
with n = m 5 (Mo = < V2/2). The critical point 

(X, Y) = (2V2/3, 1/3) corresponding to the power-law solu- 
tion (3.27) is the attractor. As for initial values of ((4/ms, 
((p/ms) = t x d(<j>/rn$)/dt, p r /m\) at m$t — 1.0, we set 
(1) (5.0 xlO 1 , 1.0 , 1.0 xlO" 2 ), (2) (5.0 xlO 3 , 1.0 , 1.0 
xlO -2 ), (3) (5.0 xlO 5 , 1.0 , 1.0 xlO -2 ), (4) (5.0 xlO 1 , 0.0, 
1.0 ), (5) (5.0 xlO 5 , -1.0 , 1.0 xlO -2 ), (6) (5.0 xlO 3 , -1.0 , 
1.0 xlO -2 ), (7) (5.0 xlO 1 , -1.0 , 1.0 xlO -2 ), (8) (5.0 xlO -1 , 
0.0, 1.0 xlO" 2 ). 

In summary, we find that for (fi/m^) 6 > 4, the power 
law solution corresponding to the critical point A is an 
attractor, while for (n/m^f < 4 the scaling solution cor- 



1. the kinetic dominant solution 

First, we show that the kinetic dominant solution 
( p. 12 ) is obtained as an attractor from the various ini- 
tial conditions (see FigJ^). The kinetic dominant solu- 
tion is described by TL = (\j\2)V 2 , and the solution of 
(fu|) is given by TL = 1/6, V = V2. (H and V are 
defined by Eq.(<h4) and Eq.(16), respectively.) Even 
though we show only the case for the model V(4>) = 
H 7 cf>~ 3 , the similar results are obtained for the model with 
V((j)) = ^ 4+Q -a , a > 2, and V((f>) = ^ 4 exp[-A^-]. 
Fig. H shows that the kinetic dominance realizes at first, 
then the solution eventually approaches the critical point 
{H = 1/6, V = \/2), 



(l) V — 



(2) 



(3)--> 



(4) 



^ ':/ : 
-. // ; 
>' - 

critical point : 
"(7) 1 2 '3' 4 

V 

FIG. 5. Global stability of the solutions for the model 
with V(<f)) = /i 7 -3 . This shows that the kinetic domi- 
nance (on the hyperbola TL = (1/12)7 :>2 ) realizes at first, 
then the solution eventually approaches the critical point 
( TL = 1/6, V — \/2), which corresponds to the solution 
( |3.12[ ). We set p = m$. As for initial conditions ((j>/m$, 
(cp/ms)' = t x d(cf) / ms) / dt) at m^t = 1.0, we choose the fol- 
lowing seven sets of initial data; (1) (4.0 xlO -1 , 0.0), (2) (5.5 
xlO -1 , -3.5), (3) (5.0 xlO -1 , 0.0), (4) (1.0, -2.9), (5) (1.0 , 
0.0), (6) (2.0 xlO 1 , l.OxlO" 1 ), (7) (1.0 xlO 1 , 0.0). 



2. the inflationary solution 



Next, we show the inflationary solution (3.4) real- 
izes from the various initial conditions numerically (see 
Fig.^). The slow-rolling condition gives TL = (1/6) V, and 
the solution (|1) is described by TL = 1/12, V = 1/2. 
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(H and V are defined by Eq.(f4.44|) and Eq.(f4.45|), re- 
spectively.) From Fig. we easily find the inflation- 
ary solution for a wide range of initial data as expected. 
Even though we show only the case for the model with 
V(<fr) = fj, the similar results are obtained for the 
model with V(<j)) = p A+a (j)~ a , a<2. 
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FIG. 6. Global stability of the solutions for the model 
with V(<j>) = /i 5 <^> _1 . This shows that the slow-rolling 
(Tt — (1/6)V) realizes first, then the solution eventually ap- 
proaches the critical poin t (V = 1/2, 7i = 1/12), which cor- 
responds to the solution (3.4). We set p — ms. As for initial 
conditions ((^/ms, (0/ms) = tx d(cj> /ms) / dt) at m^t = 1.0, 
we choose the following seven sets of initial data; 
(1) (1.0, 0.0), (2) (1.0, 2.0), (3) (1.0, 3.0), 
(4) (1.5, 3.0), (5) (3.0, 3.0), (6) (l.OxlO 1 , 3.0 ), (7) (1.0 xlO 2 , 
0.0). 



C. radiation dominated stage 



1. the power-law solution (a 7^ 2) 



Here, we show that the power-law solution (3.27) and 
(3.28) (a 2) realizes from the distant point in the pa- 
rameter space numerically (see Fig.^). The power-law 
solution is given by [J-.V) = (1, 0). Because this system 
can be written as a plane-autonomous system, and there 
are no other stable fixed point, we only show one exam- 
ple. Even though we show only the case for the model 
with V(4>) = n 7 ct)~ 3 , the similar results are obtained for 
the model with V(<f>) = 2 < a < 6. 
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FIG. 7. The attractor behavior of the power-law solu- 
tion. We set p = 7715. As for initial condition, we start from 
mst — 1.0. The initial values of (T,V) are 
(4.6 x 10 _1 , — 1.8 x 10 _1 ), which corresponds to the initial 
values of (<^/m 5 , (cf>/ms)' = t X d(<j> / m$) / dt) are (1.0, 0.0). 
Furthermore, p^/mf = 1.0 and p r /mf = 1.0. 



2. the kinetic dominant solution 



Finally, we show the result for the models with V{4>) = 
p, 4+a cj)- a , with a > 6 and V($ = ji 4 e _A ^ m5 , in the 
radiation dominant stage (Fig.|8|). Because this kinetic 
dominant solution does not include the characteristic 
mass scale, we cannot introduce the fixed critical point 
like the previous solutions. Therefore, we show the ra- 
tio of the potential energy of the scalar field to the to- 
tal energy of the scalar field. From the potential dom- 



inant initial condition [pi /p, 



(P) lAtot) 



1 initially) , the 



kinetic term eventually dominates the potential term. 

I P$ 0t ^ ~ * 0). Even though we show only the case 
for the model V(4>) = /i 11 ^ -7 , the similar results are ob- 



tained for the other models V{4>) = p 
and V{(j>) = ^ 4 e"" A< ^ m5 . 
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FIG. 8. The ratio of the potential energy of the scalar field 
to the total energy of the scalar field. This shows that even 
though the potential energy dominates the kinetic energy ini- 
tially, the kinetic term dominance eventually realizes. We set 
(j, = 7715. As for initial condition, we start from m^t = 1.0. 
The initial values of ((p/rns, (0/ms)' = t X d(<f> / 'ms) / dt) 
are (1.2,0.0). Furthermore, p^ ot) /mj = 2.8 x 10 _1 and 
Pr/mf = 1.0 

VI. SUMMARY AND DISCUSSION 

In this paper we have studied the dynamics of a scalar 
field with the typical potential in a brane world scenario. 
We have adopted the second Randall-Sundrum brane sce- 
nario and assumed that the scalar field is confined on our 
4-dimensional spacetime. In this model, because of the 
quadratic term of the energy density, which is from the 
consequence of our brane embedded in the extra dimen- 
sion, the dynamics of the scalar field is modified from a 
conventional cosmology in the early stage of the universe. 

For the power-law potential model, which was dis- 
cussed in p2fl , there exist an inflationary solution in 
the quadratic term dominant stage, which satisfies slow- 
rolling condition. Similar to the conventional cosmology, 
the e-folding number more than 60 is obtained for suf- 
ficiently large initial values of <fi. However, because the 
cosmic expansion is faster than that in the conventional 
cosmology, a shorter time is required for the enough e- 
folding number. 

After the slow-rolling condition is broken, and the os- 
cillating phase starts, the dynamics of the scalar field de- 
pends on the power of the scalar field potential. For the 
model, V(4>) = ^m 2 <j> 2 , the energy density of the scalar 
field decreases as the pressureless perfect fluid and more 
slowly than that in the conventional case. In the pre- 
heating stage, this changes the time dependence of the 
variable appeared in the Mathieu equation and leads to 
an efficient particle production via non-perturbative de- 
cay of inflaton even for the smaller coupling constant. For 
the model, V{(f>) — |A</> 4 , the energy density of the scalar 
field decreases as the radiation fluid. In this case, the so- 
lution of the scalar field evolution is represented in terms 
of the elliptic function and the form of the Lame equa- 
tion in the preheating stage is not changed. Although 
this is the same as the conventional case, because of the 
slower growth of the scale factor, the particle production 
becomes a little efficient |l8| ]. 

For the models with an inverse-power-law potential 
model (V — /i 4+Q 0~ Q ) or an exponential potential model 
V = /i 4 e~ A( ^/ m:s , we find much difference from conven- 
tional cosmology. We find analytic solutions in several 
circumstances, and analyze those stabilities. According 
to the solutions and those analysis, the aymptotic behav- 
ior of the universe is sammarized as follows: 

For V = [i 4+a <f)- a , with a < 2 and with a = 2 
(fi > \f4Qm§), an inflationary solution is obtained in the 



quadratic term dominant stage. Since this is an attrac- 
tor, such a solution is eventually obtained even though 
from the radiation-dominant initial condition. The infla- 
tionary solutions in these models are milder than that in 
the conventional cosmology. Because this potential has 
no minima, in order to recover the Big-Bang cosmology, 
some unknown reheating mechanisms is required . The 
gravitational particle production could provide its reso- 
lution |§. 

For V = pi 6 4>~ 2 with \/i ms < p < m 5 , we 

find one critical solution, i.e. a ~ t p , with p = [1 + 
1/8(/i/to 5 ) 6 ]/6. The power p is in the range of 1/4 < p < 
1. This solution is the attractor for p, > \f\ ms, which 
gives = 1 and 1/4 < p < 1. While for \i < \/i m 5 , 
another critical point with p = 1/4 appears, which turns 
to be the attractor. For this scaling solution, depends 
on (j, (n^ = t>(-^) 3 ), and the small value of p, (p < 
Zfi ms), the energy density of radiation in this attractor 
universe exceeds that of the scalar field. 

For V — p i+a <fi~ a , with 2 < a < 6, which we dis- 
cussed in the previous paper for the quintessence sce- 
nario [ p^|p^ ], the attractor universe is radiation domi- 
nant and the power-law solution is the attractor for the 
scalar field. Our work concerned to quintessence for this 
model is completed in this paper to show the stability of 
the solutions. (For the naturalness of the quintessence 
scenario and the constraints to this model, sec . ) . 

For V = p, 4+a <p- a , with 6 < a and V = e"^^ 7 ™ 5 , 
because the shape of the potential is too steep, the kinetic 
term of the scalar field dominates the potential term both 
in the scalar field dominant universe and the radiation 
dominant universe asymptotically. 

In this paper, we have assumed that a scalar field is 
confined in the brane. However, most models based on 
a string theory require a dilaton field in the bulk, which 
may gives very important effects on the behavior of grav- 
ity and a scalar field on the brane |l^]. The analysis in 
such models is left to the future work. 
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